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ABSTRACT 

The  results  of  velocity  measurements  performed  on  specie 
mens  of  polycrystalline  forsterite  with  6%  porosity,  a  single 
crystal  spinel,  and  on  obsidian  are  reported.  The  velocities 
of  compress ional  and  shear  waves  were  determined  at  both  ambi- 

G- 

enc  (25°C  and  1  bar)  conditions  and  as  a  function  of  pressure 
on  all  three  specimens.  Velocity  as  a  functior  f  tempera¬ 
ture  was  measured  to  150°C  at  1  atm  on  ‘he  obsidian  specimen. 

The  results  are  given  by  the  following  equations  for  the  two 
isotropic  materials  (forsterite  and  obsidian) ;  P  i3  in  kbars 
and  T  is  in  °C: 

Forsterite: 

vp  a  7.586  +  1.03  ±  1.03  ±  5%  x  lo"2  P 

Vg  =  4.359  +  2.45  ±  5%  x  10”3  P 

Obsid*  an: 

vp  «  5.713  -  1.87  x  10“2  P 

v s  =  3.526  -  1.43  x  10”2  P 

The  results  for  the  spinel  are  given  in  terms  of  the  single 
crystal  elastic  constants: 

C11  c12  c44 

1  bar  2985.7  1537.2  1575.8  (kbar) 

2  kbar  2995.5  1545.0  1577.5  (kbar) 
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The  isotropic  bulk  modulus  B,  Poisson's  ratio  o,  and 
the  pressure  derivatives  of  the  bulk  modulus  (dB/dP)  com¬ 
puted  from  the  above  are: 


Bs 

dBg/dP 

0 

Forsterite 

973.6 

kbar 

4.8 

0.254 

Spinel 

2019.9 

kbar 

4.18 

0.260 

Obsidian 

378.6 

kbar 

-0.87 

0.192 

The  derivatives  of  the  velocities  and  of  the  bulk  modulus  for 
obsidian  were  found  to  have  an  opposite  sign  from  those  values 
found  for  the  other  materials. 
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I.  INTRODUCTION 

The  work  performed  in  the  period  December  15,  1965,  to 
December  14,  1966,  is  presented  *n  this  report.  The  studies 
carried  out  involved  the  measurement  of  the  elastic  properties 
nd  their  variation  with  pressure.  Materials  included  in  this 
report  are  a  polycrystalline  forsterite,  single  crystal  spinel, 
and  obsidian  (glass  of  volcanic  origin) .  The  specimen  of 
forsterite  is  oi  a  rather  poor  quality,  having  about  6%  por¬ 
osity.  The  importance  of  forsterite  to  the  studies  concerning 
the  earth's  mantle  is  such  that  data  upon  a  poor  specimen  are 
useful. 

The  program  differs  from  the  usual  experiments  on  the 
elastic  properties  of  rock  in  that  we  perform  our  measurements 
upon  specimens  of  sufficient  quality  in  order  to  obtain  as 
great  an  accuracy  as  possible.  From  these  measurements  we 
have  been  able  to  define  the  pressure  derivatives  of  the  elastic 
properties.  These  derivatives  are  used  to  extrapolate  our  data, 
taken  at  very  low  pressures,  into  the  regime  of  very  high  pres¬ 
sures. 

II.  PUBLICATIONS 

During  this  year,  the  following  papers,  which  resulted 
from  our  work,  have  appeared  or  were  submitted  for  publication: 

1.  "The  Pressure  Derivatives  of  the  Sound  Velocities  of 
Polycrystalline  Alumina,"  E.  Schreiber  and  0.  L.  Anderson, 

J.  Am.  Ceram.  Soc.,  49(4],  184-190,  1966. 
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2.  "Seismic  Parameter  $:  Computation  at  Very  High 
Pressure  from  Laboratory  Data,"  0.  L.  Anderson,  Bull.  Seismol. 
Soc.  Am.,  56(3],  725-731,  1966. 

3.  "Temperature  Dependence  of  the  Velocity  Derivatives 
of  Periclase,"  E.  Schreiber  and  0.  L.  Anderson,  J.  Geophys. 

Res.,  71(12),  3007-3012,  1966. 

4.  "A  Variable  Impedance  Matching  Transformer,"  Paul 
Mattaboni  and  E.  Schreiber,  Rev.  Sci.  Instr.,  32(11) »  1625- 
1626,  1966. 

5.  "Pressure  Derivatives  of  the  Sound  Velocities  of 
Polycrystalline  Forsterite,  with  6%  Porosity,"  E.  Schreiber 
and  0.  L.  Anderson,  J.  Gecphys.  Res.,  in  press. 

6.  "Elastic  Moduli  of  Synthetic  Single  Crystal  Spinel 
at  25°C  and  to  2  Kbar,"  E.  Schreiber,  submitted  for  publica¬ 
tion. 

III.  MEASUREMENTS  AND  RESULTS 

A.  Polycrystalline  Forsterite  and  Single  Crystal  Spinel 

The  results  of  the  studies  upon  a  specimen  of  polycrys¬ 
talline  forsterite  with  6%  porosity  and  a  specimen  of  a  single 
crystal  spinel  are  reported  in  the  attached  preprints.  Appendices 
A  and  B. 

B.  Obsidian 

Velocity  measurements  have  been  performed  upon  a  sample 
of  obsidian  up  to  pressures  of  3  kbar  at  25°C,  and  the  pres¬ 
sure  derivatives  of  the  velocities  and  bulk  modulus  have  been 
determined.  The  results  of  these  measurements  indicate  that 
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the  behavior  of  the  elastic  properties  of  obsidian  differs  from 
that  which  we  have  found  for  polycrystalline  materials. 

The  specimen  which  we  used  came  from  Oregon,  and  was  pre¬ 
pared  in  the  usual  way  in  the  form  of  a  cube  about  1  cm  on 
edge.  The  chemical  analysis  of  this  sample  is  given  in  Table  I. 
The  density,  determined  by  Archimedes1  method,  was  found  to  be 
2.357  gms/cm1. 

The  velocities  were  determined  using  the  method  of  pulse 

A 

superposition,1  described  in  our  previous  reports.  The  de¬ 
pendence  of  velocity  upon  pressure  is  calculated  from  the  mea¬ 
sured  relative  change  in  frequency  (f/f0)  and  the  relative 
change  in  length  ( l/lQ )  using  equation  (1) 

v  *  vQ  <f/fc)  (*/i0)  (1) 

The  ratio  (i/i0)  is  calculated  from  the  frequency  data  using 

the  equations  given  by  Cook**  and  described  in  our  previous  re- 
o 

ports. 

The  adiabatic  bulk  modulus  (B)  and  its  pressure  deriva¬ 
tive  (dB/dP)  and  Poisson's  ratio  (o)  are  computed  from  the 
measured  velocities  using  the  following  relations: 
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TABLE  I.  Chemical  analysis  of  obsidian 


Constituent  Weight  % 


Si02 

76.19 

Al2°3 

12.61 

Fe2°3 

0.02 

FeO 

0.64 

MgO 

0.01 

CaO 

0.72 

Na20 

4.84 

k2° 

4.26 

h2o(+) 

0.85 

H20(-) 

0.21 

Ti02 

0.17 

MnO 

0.03 

p205 

0.01 

Total  100.56 

8. 

Equation  (3)  requires  the  volume  expansion#  av;  and  the 
Griineisen  constant#  y.  The  value  of  10 . 2xlO”6°C’'1  for  av  was 
obtained  from  a  measurement  of  the  linear  expansivity.  A 
value  of  y  *  0.15  was  obtained  by  first  estimating  the  Debye 
temperature  from  the  mean  sound  velocity. A‘  The  heat  capac¬ 
ity  at  constant  volume  was  obtained  from  the  Debye  tempera¬ 
ture^  and  the  chemical  composition;  the  Griineisen  constant 
was  ther  calculated  by  an  iteration  using  the  equations 


(1  +  Tavy) 

aBs 

(5) 

C  0 

v 

(6) 

The  results  of  the  measurements  at  the  ambient  conditions 
of  25°C  and  1  bar  are  given  in  Table  II. 

TABLE  II.  Elastic  properties  of  obsidian 
at  1  atm  and  1  bar 


Property 

Symbol 

Value 

Longitudinal  Velocity 

vp 

5.713  km/sec 

Shear  Velocity 

vs 

3.526  km/sec 

Bulk  Modulus 

B 

378.6  kbar 

Poisson's  Ratio 

0 

0.197 

Specimen  density  =  2.357 

3 

gms/cm 
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The  variation  of  the  frequency  with  pressure  is  shown 
in  Figures  1  and  2.  The  variation  of  velocity  was  computed 
from  these  data,  and  can  be  expressed  by  means  of  the  re¬ 
lations  (P  is  in  kbars) 

Vp  *  5.713  -  1.87xl0"2  ?  (7) 

Vg  =*  3.526  -  1. 43xl0~2  P  (8) 

The  value  of  the  pressure  derivative  of  the  bulk  mod¬ 
ulus  calculated  using  equation  (3)  is  dB/dP  «  -  0.87  kbars. 

IV.  DISCUSSION 

The  results  on  obsidian  are  of  interest  because  the  vel¬ 
ocities  are  decreasing  with  increasing  pressure.  This  is  in 

agreement  with  the  measurement  of  Birch5  and  Birch  and  Ban- 
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croft  upon  a  sample  of  obsidian  from  California.  This  re¬ 
versal  in  velocity  benavior  also  applies  to  the  variation 
of  the  bulk  modulus  with  pressure,  which  is  not  only  much 
smaller  than  the  value  of  4  which  we  have  found  for  crystal¬ 
line  solids,  but  is  also  negative.  A  negative  value  of 
dB/dP  is  in  agreement  with  the  findings  of  Bridgman2  for 
silica-rich  glasses. 

V.  FUTURE  WORK 

We  plan  to  explore  the  effects  of  temperature  upon  the 
elastic  properties  of  the  spinel  single  crystal  and  upon  the 
obsidian  specimens.  We  have  also  arranged  for  the  hot-pressing 
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of  dense  polycrystalline  forstente,  spinel,  and  enstatite. 
Future  work  on  these  specimens  is  contingent  upon  the  suc¬ 
cess  of  the  hot-pressing  effort.  Although  we  have  per¬ 
formed  measurements  upon  samples  of  two  of  these  materials, 
we  believe  it  useful  and  necessary  to  carry  out  further  mea 
surements  upon  dense  polycrystalline  specimens,  Particularly 
upon  the  forsterite,  as  our  present  specimen  contained  6% 

porosity. 

VI.  FINANCIAL  STATEMENT 
The  financial  statement  at  the  present  time  is: 

Estimated  expenditures  and  commitments  to  date  $223,503.00 
Estimated  funds  to  completion  44,717.00 

Estimated  date  of  completion  14  December  19o7 
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Figure  1 


Figure  2 


FIGURE  CAPTIONS 


Frequency  versus  pressure: 
obsidian. 

Frequency  versus  pressure: 
in  obsidian. 
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APPENDIX  A 


Pressure  Derivatives  of  the  Sound  Velocities 
of  Polycrystalline  Forsterite,  with  6%  Porosity^ 

Edward  Schreiber  and  Orson  L.  Anderson 
Lamont  Geological  Observatory  of  Columbia  University 

Palisades,  New  York 

The  determination  of  the  sound  velocities  of  polycrystallin 
forsterite  under  pressure  has  been  an  important  objective  of  our 
laboratory  at  Lamont.  We  have  obtained  numerous  polycrystalline 
samples  prepared  in  several  laboratories  and,  until  recently, 
have  not  been  able  to  pass  ultrasound  through  them.  We  finally 
obtained  a  sample  of  forsterite  which  yielded  some  weak  echoes 
at  20  mHz,  a  barely  marginal  condition  for  measuring  the  elastic 
constants.  We  were  able  to  determine  the  pressure  derivatives 
of  the  sound  velocities  on  this  sample.  The  results  should  be 
of  interest  to  the  scientific  community  concerned  with  the 
physics  of  the  earth's  interior  even  though  the  results  are  on 
a  porous  sample,  since  they  represent  the  only  such  data  exist¬ 
ing  for  forsterite.  Indeed,  it  may  prove  quite  difficult  to 
obtain  a  sample  which  is  both  ideal  and  measurable  with  precise 
techniques. 

Synthetic  polycrystalline  'forsterite'  was  obtained  in 
the  form  of  a  rod  about  3/8"  in  diameter  and  8'  long. 

^■Lamont  Geological  Observatory  Contribution  No.  00' 
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The  specimen  was  supplied  by  the  Silk  City  Ceramic  Company  of 
Hawthorne,  Mew  Jersey.  The  density  of  this  material  was  de¬ 
termined  to  be  3.021  or  94.0%  of  theoretical  density.  The 
specimen  consists  essentially  of  forsterite  (95%)  as  determined 
by  X-ray  diffraction.  A  second  phase  present  is  probably  barium 
aluminum  silicate,  but  positive  identification  was  not  possible. 
The  average  grain  size  of  the  specimen  is  about  12  microns.  In 
a  thin  section,  the  grains  appear  to  be  equidimensional  witn 
no  apparent  preferred  orientation. 

Examination  of  the  polished  surface  as  well  as  thin  sec¬ 
tion  indicated  a  porosity  of  5  i  II  consistent  with  the  mea¬ 
sured  deviation  of  the  density  from  the  theoretical  value,  the 
pore  structure  consisting  entirely  of  isolated  pores,  which  was 
further  verified  during  the  density  measurement..  No  gain  in 
weight  of  the  specimen  was  observed  after  it  had  been  in  bo-'xing 
water  for  several  hours.  Mcst  of  the  pores  (>90%)  are  equidimen¬ 
sional  and  about  16-32  microns  in  diameter.  Linear  pores  with 
moderate  aspect  ratios  were  observed.  (The  aspect  ratio  is  a 
measure  of  pore  elongation  and  is  defined  as  the  ratio  of  the  pore 
width  to  pore  length.  The  aspect  ratio  of  a  section  through  a 
sphere  would  be  unity,  and  through  a  fine  cleavage  crack,  near 
zero.)  The  aspect  ratio  of  the  pores  found  in  this  specimen  was 
between  0.2  and  0.1. 

According  to  Walsh  [1965]  and  Brace  [1965] ,  the  pressure 
required  to  close  an  elliptical  cavity  is  given  by 


P  =  EO 


(1) 
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where  P  is  the  pressure  required  to  close  the  cavity,  E  is 
the  Young’s  modulus  of  the  material  surrounding  the  cavity, 
and  0  is  the  aspect  ratio  of  the  cavity.  The  value  of  the 
Yeung's  modulus  of  olivine,  based  on  the  single  crystal  mea¬ 
surements  of  Verma  [1960] ,  is  2013  kb.  If  we  take  the  as¬ 
pect  ratio  of  the  elongated  pores  in  the  specimen  to  be  0.1, 
and  if  we  assume  the  pores  to  be  ellipsoids,  201  kb  are  re¬ 
quired  to  close  the  pores.  This  is  a  very  high  pressure,  in¬ 
deed,  and  it  would  appear  that  the  aspect  ratio  of  the  pores 
observed  in  the  specimen  is  too  large  for  pore  closure  to  af¬ 
fect  the  measurement. 

The  velocities  of  both  the  longitudinal  and  shtar  wave 
disturbances  were  measured  in  a  piece  of  the  specimen  1.0444 
cm  in  length  at  25°C.  The  surfaces  of  the  specimen  were  pol¬ 
ished  flat  to  0.1X  of  sodium  light  and  parallel  to  lxlO-4  cm/ cm. 
The  method  used  to  measure  the  velocities  and  their  derivatives 
with  pressure  was  pulse  superposition  [Schreiber  and  Andersen, 
1966 ;  McSkimin,  1961] .  A  carrier  frequency  of  10  mHz  was  used 
due  to  the  high  acoustic  loss  at  20  mHz.  This  is  the  first 
time  that  measurements  with  these  precise  techniques  have  been 
successfully  obtained  using  so  low  a  currier  frequency.  The 
pressu  ‘-generating  system  has  been  previously  described 
[Anderson  and  Schreiber,  1965] . 

The  longitudinal  and  shear  velocities  of  this  specimen  of 
forsterite,  at  1  atm  and  25°C  (ambient  pressure  and  tempera- 
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ture) ,  were  found  to  be 

vp  =  7.587  *  0.25%;  vs  =  4.359  *  0.25%  km/sec  (2) 

The  adiabatic  moduli  and  Poisson's  ratio  were  calcu¬ 
lated  using  these  velocities  and  the  measured  bulk  density 
of  3.021  gm/cm^  (see  Table  1). 

Verma  [I960]  determined  the  elastic  constants  for  an 
olivine  (density  *  3.324).  The  Voigt-Reuss-Hill  averaging 
scheme  applied  to  his  data  yields  a  bulk  modulus  of  1313  kilo- 
bars.  Using  the  method  of  Walsh  et  al.  [1965]  to  adjust  for 
the  density  difference  between  specimens,  the  two  values  of 
the  bulk  modulus  agree  to  within  10%.  The  difference  is  no 
doubt  due  to  compositional  variation  between  the  specimens. 

Typical  data  for  the  variation  of  frequency  with  pres¬ 
sure  are  shown  in  Figures  1  and  2.  As  a  result  of  the  scatter, 
the  slopes  are  known  to  about  -  5%. 

Pressure  variations  of  velocity  and  of  the  moduli  were 
computed  from  the  data  as  previously  described  [Schreiber  and 
Anderson,  1966] .  The  variations  of  the  longitudinal  and  shear 
mo'-.e  velocities  with  pressure  are  given  by  equation  (3)  ,  where 
P  is  in  kilobars. 

vp  =  7.586  +  1.03  ±  5%  x  10“2  P 

vg  «  4.359  +  2.45  ±  5%  x  10“2  P 


(3) 
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These  data  were  used  to  evaluate  the  pressure  derivatives 
of  both  the  bulk  modulus  and  of  Poisson's  ratio.  These  are 
(dBs/dP)  *»  4.8  t  0.5  and  (do/dT)  «  7.4  -  0.8  x  10“^.  The  error 
in  these  derivatives  arises  principally  from  the  error  in  slope 
of  the  velocity-pressure  data. 

The  value  of  the  Gruneisen  constant  y  *  0.97,  for  this 
sample,  was  calculated  with  the  expression 


Y 


av^s 

PC„ 


(4) 


and  the  isothermal  bulk  modulus  computed  to  be  966.6  kb  by 


Bs/Bt  -  (1  +  Tavy)  (5) 

6  o 

The  heat  capacity  Cp  =  7.889  x  10  dyne-cm/  C/gm  was  obtained 
from  Kubaschewski  and  Evans  [1958) ,  and  the  coefficient  of 
volume  expansion  of  this  specimen  av  -  24  x  10-^  °C*~*  from 
Soga  and  Anderson  [196*1 .  The  value  of  y  would  be  about  20% 
higher  for  a  nonporous  sample  due  to  the  larger  value  of  Bs 
at  zero  porosity. 

The  results  of  the  present  study  are  summarized  for  poly¬ 
crystalline  forsterite  with  6%  porosity  in  Table  1. 

We  are  aware  of  the  strong  anisotropy  exhibited  by  olivine 
in  rock  samples.  We,  therefore,  checked  our  specimen  for  indi¬ 
cations  of  anisotropy  in  two  ways.  First,  a  thin  section  was 
prepared  and  examined  under  crossed  Nicol  prisms  with  a  polar¬ 
izing  microscope.  Second,  we  measured  the  shear  sound  veloc- 


6. 


ities,  changing  the  orientation  of  the  vibration  direction. 

In  both  instances,  the  specimen  was  isotropic  in  its  behavior. 

The  effect  of  pore  closure  on  the  velocity  was  not  ex¬ 
pected  to  present  any  problem  in  this  case  due  to  the  pore 

o 

geometry.  The  fact  that  the  frequency  pressure  dependence 
is  linear  is  evidence  in  support  of  this  view.  When  cracks 
are  closing/  the  velocity  shows  a  very  large  nonlinear  rela¬ 
tion  to  pressure  (for  example/  see  Birch  (1961]).  We  believe 
that  the  effect  of  porosity  on  the  pressure  derivative  is 
not  large/  and  that  the  values  for  the  pore-free  material 
will  fall  within  the  experimental  limits  reported  here. 
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TABLE  1.  Properties  of  Polycrystalline 
Forsterite  with  6%  Porosity 

Property  Value  Units 


Density 

P 

3.021 

gm/cm 

Long.  Velocity 

VP 

7.586 

km/sec 

Shear  Velocity 

vs 

4.359 

km/ sec 

Bulk  Modulus 

Bs 

973.6 

kbar 

Poisson's  Ratio 

0 

0.254 

— 

Pressure  Deriv. Long .Velocity 

Ovp/3P)T  1 . 03xl0”2 

km/sec/kb 

Pressure  Deriv. Shear  Velocity 

Ovs/3P)t  2.45x10“2 

km/sec/kb 

Pressure  Deriv.  Bulk  Modulus 

OBs/3P)t  4.8 

— 

Pressure  Deriv. Poisson' s  Ratio 

<3c 

[/3P>X  7.4xl0"4 

kb~^ 
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FIGURE  CAPTION* 

FIG.  1  Frequency  versus  pressure:  Longitudinal  mode 
in  i^SiO^, . 

^IG.  2  Frequency  versus  pressure:  Shear  mode  in 
Mg2Si04. 
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APPENDIX  E 


Elastic  Moduli  of  Single  Crystal  Spinel  at  25°C  and  to  2  Kbar 

Edv?ard  Schreibe** 

Lamont  Geological  Observatory  of  Columbia  University^ 

Palisades,  New  York 

ABSTRACT 

The  elastic  properties  of  a  synthetic  single  crystal  spinel 
with  a  composition  of  ”gO*2.6  A1203  have  been  determined  at  25°C 
and  up  to  2  kbar  by  measuring  propagational  velccitites  of  ultra¬ 
sonic  waves  using  the  method  of  pulse  superposition.  The  iso¬ 
tropic  moduli  and  Poisson's  ratio  were  computed  from  the  single 
crystal  values  using  the  Voigt-Reuss-IIill  approximation. 

Values  of  the  elastic  moduli  at  1  bar  and  2  kbar  and  at 
25°C  are 


C11  c12 
1  bar  2985.7  1537.2 


1575.8  (kbar) 


2  kbar  2995.5  1545.0 


1577. c  (kbar) 


The  bulk  modulus  was  computed  to  be  2019.9  kbar,  and  che 
pressure  derivative  of  the  bulk  modulus  was  found  to  be  4.18. 
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i .  introduction 

The  purpose  of  this  communication  is  to  present  elastic 
moduli  measurements  made  upon  a  specimen  of  synthetic  single 
crystal  spinel  at  25°C  from  1  bar  to  2  kbar,  and  the  applica¬ 
tion  of  these  data  to  calculate  the  pressure  derivatives  of 

l 

the  elastic  moduli  and  associated  bulk  properties.  Since 

2  3 

elastic  moduli  data  for  both  HgO  and  AljO^  are  available, 
similar  data  for  spinel  are  of  interest.  The  spinel  for  which 
data  are  reported  here  has  the  formula  I!g0*2.61  A^O-j.  The 
value  of  the  pressure  derivative  of  the  bulk  modulus,  dB/dP, 
is  a  significant  parameter  in  evaluating  equations  of  state. 

It  has  been  shown  that  the  value  of  dB/dP  measured  at  low  pres- 
sure  allows  one  to  estimate  volume  at  very  high  pressure. 

II.  PEASUREPENT 

A  single  crystal  spinel  boule  grown  by  flame  fusion  was  ob¬ 
tained  from  the  Commercial  Pineral  Corporation  of  New  York  City. 

The  density  of  a  specimen  cut  from  this  sample  was  measured  by  a 
volume  displacement  method,  talcing  into  account  all  buoyancy  and 
surface  tension  corrections,  and  found  to  be  3.6193  -  2xl0“4 

gms/cm  .  This  value  is  1*  greater  than  the  X-ray  density  calcu- 

5 

lated  using  the  data  of  the  N.B.S.  for  spinel  (r-gAljO^)  prepared 
by  sintering  a  stoichiometric  mixture  of  the  pure  component  oxides. 
A  chemical  analysis  of  this  specimen  indicated  a  composition  of 
PgO‘2.61  A^Oj.  Since  it  is  possible  that  Al203  is  soluble  in 
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the  spinel  'structure,  thereby  raising  the  density  of  the  spinel 

6  7 

phase  (Clark  et  al.  ,  Roy  et  al.)  an  X-ray  pattern  v/as  obtained 

v/ith  our  specimen  which  indicated  that  only  a  spinel  phase  was 

present  with  a  lattice  spacing  of  7.98^.  This  spacing  corre- 

6 ' 7 

sponds  to  a  composition  '  of  IgO‘2.6  A^O-j  and  is  in  agreement 
with  the  results  obtained  from  the  chemical  analysis. 

The  specimen,  cut  into  an  oriented  cu^>e  v/ith  four  (110) 
faces  and  two  (100)  faces,  was  polished  flat  to  within  0.1A  of 
sodium  light  and  the  faces  parallel  to  1’xlCT*  in/in.  The 
crystallographic  orientations  were  determined  and  maintained 
by  X-ray  methods.  The  elastic  constants  were  determined  by  mea¬ 
suring  the  propagation  velocity  of  ultrasonic  waves  through  the 

crystal,  employing  the  interferometric  technique  of  pulse  super- 
8  9 

position.  '  This  method  has  the  concurrent  advantages  of  high 
accuracy  (0.01%)  in  measuring  the  absolute  velocities  and  great 
sensitivity  for  detecting  the  small  changes  of  velocity  with 
pressure.  This  particular  technique  of  measurement  has  been 
described  '  '  and  will  only  be  briefly  outlined. 

The  output  from  a  pulsed  oscillator  (tone  burst)  is  applied  to 
an  appropriate  quartz  transducer  which  is  bonded  directly  to  the 
specimen.  The  frequency  of  the  carrier  wave  in  the  tone  burst  is 
the  same  as  the  transducer  resonant  frequency,  and  the  repetition 
rate  of  the  tone  burst  is  varied  until  its  period  is  equal  to  the 
round  trip  lelay  time  within  the  specimen.  Under  these  circum- 
stanc  c,  the  transmitted  and  received  pulses  overlap,  i >  the  trans- 


4  . 


mitter  is  gated  in  order  to  view  the  returning  echo  pattern  on  an 

oscilloscope.  The  echoes  are  nade  to  overlap  and,  by  adjustment 

of  the  repetition  rate,  are  made  to  interfere  constructively.  At 

the  condition  of  constructive  interference,  the  following  equa- 
.  .  8  9 

tion  applies  ' 


t  *  1/P  +  <t>/2 Tff  -  n/f  (1) 

where  t  is  the  round  trip  travel  time  of  the  wave  within  the 
specimen,  F  is  the  repetition  rate  of  the  applied  signal,  f  is 
the  carrier-wave  frequency  (40  mHz) ,  4>  a  phase  shift  introduced 
at  the  transducer-bond  interface,  and  n  the  integral  number  of 
cycles  of  displacement  of  the  carrier  wave  in  overlapping  echoes. 
As  the  pressure  is  varied,  the  phase  shift  is  maintained,  constant 
by  adjusting  the  carrier  frequency  to  the  transducer  resonance 
as  it  changes  with  pressure;  thus  the  primary  measurement  is 
the  repetition  frequency  F.  The  variation  of  the  elastic  moduli 
with  pressure  is  determined  by  measuring  the  change  in  the  repe¬ 
tition  rate,  and  the  velocity  of  a  particular  vibrational  mode 
is  calculated  with  the  relation 

Vi  =  voi  <F/Foi)  U/l0)  (2) 

where  the  zero  subscripts  refer  to  the  ambient  conditions  of  25°C 

and  1  bar.  The  length  ratio  ( i/iQ )  is  calculated  from  the  fre- 

10 

quency  ratios  following  the  method  described  by  Cool..  For 
crystals  of  the  isometric  system,  the  length  ratio  if,  independent 


5. 


of  crystal  orientation. 

1 1 

As  shown  by  Hason,  the  elastic  moduli  nay  be  determined 
from  the  velocity  measurements  in  the  100  and  110  crystallographic 
directions.  In  Table  I,  these  velocities  are  listed  and  their 
relations  to  the  three  independent  elastic  moduli  are  indicated. 
Since  there  are  three  independent  elastic  moduli,  only  three  vel¬ 
ocities  need  be  measured  to  completely  define  the  system,  but 
two  extra  velocities  are  measured  for  the  purpose  of  a  cross-check. 

TABLE  I.  Relation  of  measured  velocities  to  the  elastic  moduli. 


Propagation  Direction  Velocity  and  Vibration  Formulae  for 

Direction  Velocities 


The  pressure  system  employed  has  been  described  in  Schreiber 
3 

and  Anderson  and  is  so  designed  that  pressure  and  temperature 
are  both  measured  in  situ.  ?  pressure  balance  is  employed  in 
order  that  the  precision  of  the  pressure  measurement  be  consonant 
with  the  precision  of  the  velocity  measurement.  The  vessel  tern- 


perature  is  maintained  at  25°C  by  circulating  thermostated  water 
through  a  copper  ceil  surrounding  the  vessel.  All  measurement 5: 
were  performed  within  *0.05°C  of  25°C  and  with  the  piston  of  the 
pressure  balance  floating. 

III.  RESULTS  AND  DISCUSSION 
A.  'Measurements  at  1  bar 

The  results  of  the  measurement  of  the  wave  velocities  at  the 

ambient  conditions  of  25°C  and  1  bar  are  given  in  Table  II.  Verma 

measured  the  velocities  in  a  single  crystal  spinel  and  his  results 

are  included.  A1J  High  the  composition  of  his  spinel  is  reported 

to  be  flgO'3.5  Aljt^,  there  is  close  agreement  between  the  two  sets 

of  velocities.  It  would  appear  that  substitution  of  A1  for 
2+ 

Mg  in  the  spinel  lattice  does  not  have  any  great  effect  on  the 
velocities. 


TABLE  II.  Measured  wave  velocities  in  spinel  at  1  bar. 


Designation 

Velocity  (km/sec) 

This  Work*  Verma* 

V1 

9 . 0833 

9.10 

v2 

6.5983 

6.61 

v3 

10 .  '"6 

10.30 

v4 

4.4733 

V5 

6.5978 

—  - - - - .  - - - - H - 

1  — ■  —  ■  - 

♦Densities  are  3.63  gm/cirr  for  Verma’ s  specimen  and  3.6193  gm/cm 
for  the  specimen  reported  here. 
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The  elastic  moduli  computed  from  the  measured  density  and 
velocities  are  given  in  "able  III.  The  moduli  are  computed  from 
different  combinations  of  the  wave  velocities  in  order  to  cross¬ 
check  the  measurements.  The  compliances  are  computed  from  the 
average  value  for  the  corresponding  moduli.  The  isotropic  bulk 
(B)  and  shear  (G)  moduli  are  calculated  from  the  single  crystal 
values  according  to  the  v/cll-known  schemes  of  Voigt  (for  mod  Li) 
and  Reuss  (for  compliances)  and  are  finally  reported  as  the 
averaged  Voigt-Reuss-Hill  value  following  Anderson.1 

B.  Measurements  under  Pressure 
The  change  in  repetition  frequency  with  pressure  is  shown 
in  Fig.  1 ,  where  the  frequency  ratio  is  plotted  as  a  function  of 
pressure.  The  solid  lines  represent  the  least  square  fit  to  each 
set  of  experimental  data.  The  frequency  ratios  may  be  expressed 
by  the  following  relations  (P  is  in  kbars) : 

Fl/F01  =  1-000  000  +  7.21  x  10"4  P 

F2/F02  *  1-000  000  +  2.13  x  10"4  P 

,  -4  (3) 

F3/F03  *  1-000  000  +  5.98  x  10  H  P 

F4/F04  =  1.000  000  +  2.64  x  1C“4  P 

The  results  obtained  from  the  pressure  data  are  summarized 

in  Table  IV.  The  frequency  ratios  given  by  Eqs.  (3)  were  used  to 

l  o 

determine  the  length  ratio  (following  Cook  )  and  the  density 
ratio  was  calculated  from  the  initial  density  (3.6193  gm/cm3)  and 
the  length  ratio.  The  isotropic  moduli  xrere  computed  in  the 
manner  described.  The  pressure  derivative  of  the  bulk  modulus 
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TABLE  III.  Elastic  properties  of  spinel  at  25°C  and  1  bar. 


Prope  rty 

Relation 

Value 

C11 

pvf 

2936.1  kbars 

C11 

Average 

p{vf  +  -  v§) 

2985.4 

2985.7 

c12 

P(v§  -  v|  -  v§) 

1536.9 

c12 

Average 

p(vjr  -  2v|) 

1537.6 

1537.2 

c/»  4 

pv^ 

1575.7 

c44 

Average 

pv| 

1575.8 

1575.8 

S11 

<cll+c12)/(cifcl2)  (c11+2c12) 

5.1525xlG~4  kbars"1 

S12 

"c12/<Cll”C12) (cix+2c12) 

-1.7511xl0~4  ” 

S44 

1/C44 

6.3467xl0“4 

Bv 

1/3 (c11+2c12) 

2020.0  kbars 

Gv 

1/5 (cii”ci2+3c44^ 

n  -1 

1235.0 

Br 

[3(s11+2s12)  J 

201^.5 

Gr 

[V5(8u-s12)+3/58  ^y1 

1071.7 

Bvrh 

1/2(BV+Br) 

2019.9 

Gvrh 

1/2(GV+Gr) 

1153.3 

o 


( 3B-2G) /2 ( 3B+G) 


0.2601 
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TP OLE  IV. 

Summary 

of 

results  for 

sninel  at  25°C. 

Length  ratic 

at 

?  kb  (i/lQ) 

=  0.999664 

Density  ratio 

at  2  kb  (p0/p) 

*  1.001007 

1  bar 

2  kb 

3 

Density  (gm/cm  ) 

3.6103 

3.6229 

vi 

9.C  33 

9.093^ 

V2 

6.5933 

6.5989 

Velocity  (km/sec) 

V3 

10.206 

10.305 

V4 

4.4733 

4.4742 

V5 

6.597C 

6.5904 

C11 

2905.7 

2995 . 5 

Ave. Moduli  (kbar) 

CX2 

1537.2 

1545.0 

c44 

1575.0 

1577.5 

S11 

5.1525xl0”4 

5.1438xl0”4 

Compliances 

(kbar”1) 

S12 

s.;4 

-1.7511xl0"4 

6 .  3467xl0~4 

-I.7503xl0”4 

6.3391xl0”4 

Pv 

2020.0 

202C.5 

Bulk  r.odulus 

Br 

2019.8 

20?8.5 

(kbars) 

i 

Bvrh 

2019.9 

2028.5 

Gv 

1235.0 

1236.6 

Shear  i’odulus 

Gr 

1C71.7 

1073.1 

(koars) 

Gvrh 

1153.3 

1154.8 

0.260  0.260 
Pressure  derivative  bulk  rot ‘ulus 


(AVAP)T  =  4.2  (d.B/dP)  „  = 


4.18 
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(dB'dP)  was  computed  ir  two  ways,  first,  as  ( AP/AP)  from  the 
tabulated  data  which  involve  the  approximate  length  charge 

computation,  and  second,  fro?  the  relation  derived  for  the 

13  14 

natural  velocities  by  Thurston  and  by  Thurston  and  Brugger. 

Their  relation  for  crystals  of  cubic  symmetry  is 


(dB/d.P)  =  ^C11(F1/F01)  -  </3(cn  ~  C12)  (F^/Fp^)  +  1/3  (1+Tay)  (4) 

2 

Anderson  has  pointed  out  that  the  latter  represents  computation 
Ly  a  slope  whereas  the  former  represents  computation  by  a  chord. 
The  aaree^er.t  between  the  tv  to  values,  *.20  and  *.lr’,  obtained  by 
the  two  methods  of  commutation  may  be  taler  as  evidence  for  the 
linearity  of  the  bulk  modulus  pressure  relationship. 

It  is  useful  to  compute  the  bull  clastic  properties  as  iso¬ 
thermal  values  for  two  reasons.  First,  it  is  the  isothermal 
value  which  is  obtained  in  experiments  involving  static  compres¬ 
sion?  second,  equations  of  state  used  to  predict  pressure  volume 
relationships  at  high  pressures  ordinarily  irvolvc  the  isothermal 
value  of  L  and.  of  (d.B/dF)  .  The  equation  for  the  convcrsi  n  is 


P.T  =  :  (1  +  Toy) 

Pr  =  1999  '  bar 

whore  the  T  subscript  denotes  the  isothermal  value,  T  is  t’^e 
absolute  temperature,  a  the  coefficient  of  vclure  expansion, 
and  Y  the  Crureiscr.  constant,  which  "’ay  be  computed  from 

e  Ha 

Y  c 

Y  = 


(5) 


(6) 


1.53 


11. 


where  Cp  is  the  heat  capacity  at  constant  pressure  and  p 

the  density.  Tor  the  pressure  derivative  of  the  bull:  modulus, 

the  relationship  is  somewhat  more  invcl  ed  and  is  given  by 

l  5 

(after  Overton  ) 


dBT/dP)T  = 


(dE/dP)T  +|Ta(BT/B^J  |{-2/aUT)  (dB/dT)p  -  2(dB/dP)T~|  + 

^Ta(BT/pJ|2  r(dB/dP)T  -  (1/a2)  (da/dT)p“-J 


(7) 


(dBT/dP)T  =  4.12 


The  values  used  in  evaluating  the  isothermal  properties  from  the 

adiabatic  at  T  *  ?98°K  are:  p  =  3.6193  g/cm3  and  B  =  2019.9  kbar 

(this  v;ork) ,  a  «*  22.  SxKf^Oc"1,  anddo^dT  5xl0‘8  (Skinner16), 

0^  =  27.8cal-deg’’1-mol~1  (Kubaschewski  ard  rvans17),  and  (dBT/dT) 

was  estimated  to  be  -0.21  kb°C  from  the  work  reported  by 
1  8 

Soga  et  al. 

4 

Anderson  has  shown  that  for  the  condition  where  the  bulk 
modulus  is  linear  with  pressure  the  Turnaghan  equation  of  state 
integrates  to  (V  is  volume) 

*n(V0/V)  -  j^l/(dP,T/dP)T^J  in  £(dBT/dP)T(P/BT)  +  lj  (8) 

which  contains  no  arbitrary  constants  and  he  has  demonstrated  the 
applicability  of  Iq.  (8)  to  numerous  compounds.  The  agreement 
between  the  two  methods  of  computing  (dB/dP) ,  shown  in  Table  IV, 
indicates  that  BT  is  linear  with  pressure  and  therefore  Eq.  (8) 
should  apply  for  spinel.  In  this  case,  the  equation  of  state 


for  spinel  is 


12  . 


V  r  -3  1  0.239 

-2.  =  |_  2.06x10  P  +  1 J  (9) 

Equation  (9)  should  hold  up  to  at  least  1000  kbars  if  there 
is  no  phase  change. 
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FIGURE  CAPTION 


Fig.  1 


Variation  of  frequency  ratio  with  pressure. 
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